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ABSTRACT 


An analytical study was conducted to determine approximate radiator transfer func- 
tions for use in SNAP- 8 ground test facilities. Theoretical equations for transient heat 
transfer and flow processes in the radiator were programmed on a digital computer. 
With this computer simulation, the frequency response of the radiator -outlet tempera- 
ture was determined. Approximate transfer functions were determined from the 
frequency-response information. The transfer function determined for outlet tempera- 
ture to inlet temperature was a single lag of third order cascaded with a dead-time func- 
tion. The time constant and dead time were inversely proportional to flow rate. The 
simplified transfer function chosen to describe the response of outlet temperature to in- 
let flow rate was a first-order lag. The first-order lag time constant was in inverse 
proportion to the flow rate. 



TRANSFER FUNCTIONS FOR RADIATOR DYNAMICS IN THE SNAP-8 SYSTEM 


by Robert W. Leko and Kent S. Jefferies 
Lewis Research Center 

SUMMARY 

In the design of practical controls for space- radiator simulation, simple mathemati- 
cal representation of radiator dynamics is needed. An analytical study, therefore, was 
conducted to determine approximate radiator transfer functions for use in SNAP- 8 ground 
test facilities. Theoretical equations for transient heat-transfer and flow processes in 
the radiator were programmed on a digital computer. With this computer simulation, 
the frequency response of the radiator-outlet temperature was determined for small per- 
turbations in inlet temperature and flow rate around steady-state values. Approximate 
transfer functions were determined from the frequency-response information. 

The transfer function determined for outlet temperature to inlet temperature was a 
single lag of third order cascaded with a dead-time function. This form agreed very 
well with the computer -generated frequency response for a range of flow rates from de- 
sign value to 10 percent of design value, provided the time constant and dead time were 
varied inversely with the flow rate. 

The frequency response generated by the computer program for the radiator outlet 
temperature to flow rate was like that of a first-order lag with a resonance effect super- 
imposed in the midfrequency range. The simplified transfer function chosen to describe 
the response of outlet temperature to inlet flow rate was a first-order lag. The first- 
order lag form of transfer function adequately covered flow rates from design to 10 per- 
cent of design value, provided the time constant was adjusted in inverse proportion to 
the flow rate. 


INTRODUCTION 

SNAP- 8 is a nuclear electrical generating system currently under development for 
space power applications (ref. 1). The power- conversion system of SNAP-8 uses the 
Rankine cycle and transfers the waste heat of this cycle to a heat-rejection loop. The 
working fluid of the heat rejection loop is sodium-potassium eutectic (NaK). The heat- 



rejection loop uses a radiator to reject the waste heat to space. For ground tests of 
systems it is usually necessary to use convectively cooled heat exchangers which have 
heat transfer characteristics substantially different from those of a flight-type radiator 
operating in space, where heat rejection is strictly radiative. Thus, in order to conduct 
meaningful studies of system dynamics in a test cell, some method of forcing the convec- 
tively cooled heat exchanger to behave like a flight radiator is required. Therefore, a 
radiator simulator is employed in ground testing of SNAP-8. 

The design and performance of a radiator simulator for a SNAP- 8 test facility at the 
Lewis Research Center is discussed in reference 2. In this reference, it was brought 
out that, for improved design of radiator simulators, more definition of the radiator dy- 
namics is required in the form of simple mathematical models that are compatible with 
the limited control equipment available in test facilities. Therefore, an analytical study 
was conducted to determine simple transfer functions that approximate the dynamics of 
the SNAP- 8 radiator. 

The analysis consisted of deriving the appropriate heat balance and flow equations 
from which a dynamic digital- computer simulation of the radiator was programmed. A 
cylindrical radiator 9. 26 meters in length and 3. 72 meters in diameter was assumed for 
this analysis although no final radiator has been selected at present for SNAP- 8. The 
computer simulation was used to generate frequency response plots of radiator outlet 
temperature to inlet temperature at various flow rates. Also, frequency response was 
obtained for outlet temperature to flow rate at various steady- state values of flow rate. 
From the frequency-response maps, the desired transfer functions were determined. 

The report includes the derivation of the equations and the computer program in addition 
to presenting the frequency response and resulting transfer functions. 


SYMBOLS 


Mathe- 

matical 

symbol 



FORTRAN Description 

symbol 

o 

cross-sectional area, m 

p 

area of fin for heat transfer by radiation, m 

p 

area of tube armor for heat transfer by radiation, m 
specific heat of tube and armor, J/((kg)(K)) 
specific heat of fin, J/((kg)(K)) 
specific heat of fluid, J/((kg)(K)) 
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Mathe- 

FORTRAN 

Description 

matical 

symbol 

symbol 


E 


internal energy, J/kg 

H 


enthalpy of fluid, J/kg 

L 


time lag, or dead time, sec 

M 


exponent for transfer function for temperature disturb- 
ances, dimensionless 

m A 


mass of tubes and armor, kg 

m F 


mass of fin, kg 

nij 


mass of fluid, kg 

N 


lump or node number, dimensionless 

Qi 


heat transferred to fin, J/sec 

q 2 


rate of energy transfer from fluid to tube by forced 
convections, 6. 28 ur (t - T)/A cg , J/((m^)(sec)) 

R 


convective heat transfer resistance, ((sec)(K))/J 

r 


tube radius, m 

T 

TDRA 

temperature of fluid, K 

t n 


temperature of fluid at N** 1 node, K 

t f 


temperature of fin, K 

T d 

= 619 

design temperature of fluid, K 

t 

TRA 

tube temperature, K 



tube temperature of N** 1 lump, K 

U 


o 

overall heat transfer coefficient, J/((m )(sec)(K)) 

V 


fluid velocity, m/sec 

W 


mass of fluid in lump, kg 

w 


fluid flow rate, kg/sec 

w d 

= 4. 54 

design point fluid flow rate, kg/sec 

AX 


lump length, m 

1/Rcw 


6. 28 rUW/pcwA cg , dimensionless 
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Mathe- 

matical 

symbol 

e 

a 


P 

9 

6 


r D 

T o 

T l> t 2’ t 3 
t 4 

a e 

D/D0 


FORTRAN Description 

symbol 


DTIME 

TIME 

C7 

FLOWD 
A159 = 170. 

TRAD 

FINEF 

DFIN14 

AMPL 

OMEGA 

STEP 

AMPLR 

RETOT = 0.02184 


emissivity 

2 4 

Stefan- Boltzmann constant, J/((m )(sec)(K )) 

Q 

density, kg/m 
time, sec 

AX/v = W/w, time for fluid to move through one lump, 
sec 

time constant, sec 

design point time constant, sec 

time constant of third-order transfer function, sec 

general time constants, sec 

time constant of first-order lag transfer function, sec 

multiplication sign 

calculation time interval, sec 

substantial derivative 

time, sec 

inlet fluid temperature, K 
flow rate per tube, kg/sec 
conductivity of armor, J/((m)(sec)(K)) 
average radiator surface temperature, K 
fin effectiveness 

derivative of the fourth root of fin effectiveness with 
respect to radiator surface temperature, K~* 

amplitude, K or kg /sec as appropriate 

frequency, rad/sec 

step change variable, K or kg/sec as appropriate 

amplitude ratio, dimensionless 

thermal resistance of tube wall, ((sec)(K))/J 
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Mathe- 

matical 

symbol 

FORTRAN 

symbol 

Description 


CON(l) = 233. 3 

reciprocal of mass of fluid per node per tube, 1/kg 


CON(2) = 0.05212 

thermal conductivity multiplied by area and divided by 
specific heat of fluid, kg/sec 


Cl = 5.103X10" 8 

emissivity of radiator surface multiplied by the Stefan- 

2 4 

Boltzmann constant, J/((m )(sec)(K )) 


HCAPTA = 35. 26 

heat capacity of tube and armor, J/K 


HCAPF =1.977 

heat capacity of fin, J/K 


AARM = 0.001918 

radiation area of armor surface for one lump of one 
tube, m^ 


AFIN = 0.015 

radiation area of one lump of two fins, two fins per 
tube, m^ 


FINWI = 0.03414 

fin width, m 


FINTH = 0. 00762 

fin thickness at root, m 


COMPUTER SIMULATION 
Derivation of Equations 

The physical description of the heat-rejection loop radiator model assumed in this 
analysis is shown in figure 1. It is a cylindrical shape 9. 26 meters in length and 3. 72 
meters in diameter. The radiator model assumed is comprised of 153 tube-and-fin 
segments which make up the radiating surface. The heat passes from the NaK fluid 
to the tube wall by forced convection and is then conducted from the tubes to the fins 
where it is radiated to space. 

In the analysis it was assumed that any one tube and fin combination behaves exactly 
like all the others so that we may focus our attention on just one for simplicity. Next, 
axial heat conduction in the tube wall has not been taken into account. The analysis is 
one dimensional. Variation in material transport properties as a function of tempera- 
ture has not been taken into account. There are no solar, planetary, or other external 
heat source inputs to the radiator assumed in this analysis because no SNAP- 8 mission 
has been defined. 
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In writing the computer program it was convenient to divide the radiator into axial 
lumps as shown in sketch (a) . 

The energy balance for a unit volume of the fluid can be stated as follows: 


Rate of accumulation 


Net rate of internal 


Net rate of energy 

of internal energy 

+ 

energy convected 

= 

transferred by 



with stream 


forced convection 


In equation form this energy balance appears as 

d(pE) , a(pvH) _ q (!) 

ae ax 2 

pc f dT [ vc dT _ 6. 28 Ur(t - T) 
d0 f dX _ A cg 


or 


= ~6. 28 Ur(T - t) (substantial derivative) 
d0 pc f A cg 


or 
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DT _ -6. 28 Ur 


( 2 ) 


T " t P° f A cs 

Assuming that the tube temperature is constant for each lump and during the time period 
from 9 - 6 to 9, the previous equation was integrated to give the following equation: 

r -,-6.28 UrW/wpCpA 

T N (0) " ' 6) = L T N-1 (0 " 6) * ^-1^ " 6) J e (3) 

Equation (3) was rearranged to give the following equation for calculating fluid tempera- 
ture: 


T N (0 ) = T (N-1)( 0 - 8) + (l - e 1/RCfW )jt( N1) (0 - 5) - T (N _ 1} (fl - ojj (4) 

The operational- form of this equation as it appears in the computer program is 

TDRA(K) = TDRA(K - 1) + FACT1 * (TRA(K - 1) - TDRA(K - 1)) (5) 


An energy balance was made over the tube and armor in the manner shown below. 


Heat transferred 


Change in internal 


Heat energy radiated 


Heat transferred 

from fluid to tube 


energy of tube 


by tube to space 

+ 

to fin 


This energy balance can be written in equation form as 

/ -1/Rc f w\ dt N 4 

< T N -t N )(l-e jc f w - m a c a — = «V N + Q 1 <6) 

An energy balance was made over the fins in the following manner: 


Heat transferred 


Change in internal 


Heat radiated 

to fin 


energy of fin 

+ 

by fin to space 


In equation form the above energy balance as applied to the fins appears as 
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( 7 ) 


dT F 

Ql — HlriCri + (JcA-pTpi 

1 F F dB F F N 


Using the definition of fin effectiveness as the ratio of heat radiated at the average fin 
temperature T p to the heat that could be radiated if the fin were at the fin root temper- 
ature t. 


aeApTj, 


fi = 


N 


aeAptjj 


or 


= fi 1 / 4 . 


N 


l N 


and- 


dTp N _ 0 l/4 dt N t dfi 1//4 dt N 
d0 dd dt N d0 

Equation (7) can then be rewritten in terms of fin root temperature (tube temperature) as 
follows: 


Ql = m F c F [0 


1/4 t dfi X / 4 \ dt N A nt 4 
+ t N — — }— + aeA F 5it N 


J 


( 8 ) 


Combining equations (6) and (8) 


/ - 1 /Rc f w\ 

< T N -t N )(l-e f ) 


dt. 


c f w - m a c a -jf = « A t l N 


+ m F c p 


^1/4^ dSJ^)^N +oeA ^4 (9) 

\ N dt N J de F N 


Rearranging equation (9) 
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( 10 ) 


/ -1/Rc f w\ ( 4 4 x 

dt N < T N " W y 1 “ 6 ] C f W “ CT£ (Vn + A F fit N/ 


d0 


m a c a +m F c Fr T L N 


1/4 . df^ 

+ ^XT 


dt 


N 


Integrating equation (10) over the time interval (0 - Ad) — 6 . 


^n( - tjq-(0 ” A0) + 


(T N " ^ 


/ -1/Rc f w\ 4 

11 - e Jc f w - ore(A t + fiA p )t N 


m A c A + m F c F( fi1 ^ 4 + ~ 


l/4> 


dt 


N 


A 9 (11) 


The operational form of equation (11) as it appears in the program is 
TRA(J) = TRA(J) + ((TDRA(J) - TRA(J)) * FACTOR/RETOT - HEAT( J))/(HCAPTA 

.+ HCAPF *' (FIN14 + TRA(J) * DFIN14)) * DTIME (12) 
Equation (12) solves for the temperature of the tube at each lump. 


Computer Techniques 

The computer program consisted of two parts: (1) a subroutine subprogram, and 
(2) a main program. 

The operational equations (eqs. (5) and (12)) which appeared previously are con- 
tained in the subroutine. A complete listing of the subroutine subprogram is included in 
the appendix. 

The main program provides for the initialization of constants and for other informa- 
tion required by the subroutine subprogram. The main program makes a series of calls 
to the subroutine. During this series of calls the subroutine then solves for the steady- 
state temperature profile in the radiator and fluid for a given set of conditions. Using 
the steady- state condition returned to the main program from the subroutine, the main 
program makes a second series of calls to the subroutine to obtain the dynamic effects. 
This is accomplished by means of a sinusoidal analysis. In the case of inlet tempera- 
ture disturbance, the perturbing function is a sine wave in inlet fluid temperature. The 
inlet temperature perturbation equation in the main program is 
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TDRA(l) = 619. 0 + AMPL * SIN(OMEGA * TIME) + STEP 


(13) 


The value of amplitude was 5. 56 K for flow rates of 10, 20, 50, and 100 percent of de- 
sign flow rate. The sinusoidal analysis is achieved by setting amplitude and omega to 
some appropriate values and STEP to zero. A step change can be simulated by setting 
amplitude and omega to zero and assigning an appropriate value to STEP. 

The analysis in the case of flow rate disturbance proceeds in an exactly analogous 
manner to that previously described for inlet temperature disturbance. Using the 
steady- state condition returned to the main program from the subroutine, the dynamic 
effects induced by a change of inlet flow rate is programmed as the equation given below: 

FLOWD = (4. 54 + AMPL * SIN(OMEGA * TIME) + STEP/153. 0 (14) 

The value of amplitude is 0. 00693, 0. 01375, 0. 0346, and 0. 0693 kilogram per second 
for flow rates of 10, 20, 50, and 100 percent of design flow rate. The values were 
chosen to correspond with the same percentage change of amplitude as in the inlet tem- 
perature perturbation studies. A complete listing of the main program is included in the 
appendix. 


FREQUENCY RESPONSE AND TRANSFER FUNCTIONS 

The dynamic response of the SNAP-8 radiator exit temperature to a change in inlet 
fluid temperature can be found in terms of its frequency response characteristics. The 
Bode plot shows amplitude ratio and phase angle plotted against frequency. The ampli- 
tude ratio against frequency portion of the plot is made on log-log paper. This enables 
adding amplitude ratios directly to build up more complex structured transfer functions. 
The phase angle against frequency part of the plot is done on semilog paper. Phase 
angles may also be added directly to build up more complex structured transfer func- 
tions. 

A parametric study was made of the radiator response of exit temperature to a 
sinusoidal change in inlet temperature. The parameter was fluid (NaK) flow rate. Data 
for plots were generated for 10, 20, 50, and 100 percent of the design point flow rate 
value. The design point flow rate is 4. 54 kilograms per second. The design point fluid 
inlet temperature is 619. 0 K. The frequency response data was generated for a constant 
inlet temperature. 

The Bode plot (fig. 2) shows the dynamic response of the SNAP-8 radiator exit tem- 
perature to changes in inlet temperature. Throughout this report, the same symbol is 
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used consistently to indicate a particular value of the flow rate parameter. A circle, 
for example, indicates design point flow of 4. 59 kilograms per second. The amplitude 
ratios are normalized; that is, the amplitude ratio at each frequency is divided by the 
steady- state value of amplitude ratio. 

One criterion of this study was to obtain the simplest form of transfer function ade- 
quately representing the radiator dynamics. This criterion arose because it was desired 
to have the simplest analog computer in a ground test simulator for controlling the radi- 
ator dynamics. A transfer function of the form shown below fits the Bode plots of fig- 
ure 2 to give the proper amplitude ratio over the frequency range studied 


'at 


out 


AT. 

in 


1 

(1 + TjSKI + t 2 S)(1 + r 3 S) 


Consistent with the criterion for simplicity, a transfer function of the form shown below 
was selected instead. 


AT in (1 + t q S) 3 

This result agrees favorably with the conclusions given in reference 2. 

The phase angle portion of the Bode plot shown in figure 2 dictated the need for a 
time lag, or dead time, element. The time lag transfer function is given below. 

„ -LS 

G = e 


The total transfer function as indicated by the Bode diagram therefore becomes 


'at 


out 


-LS 


AT in (1 + t q S) 3 


The solid lines shown in figure 2 represent the locus of points given by the simpli- 
fied transfer function. The important result was found that the time constant r Q is 
given by the relation shown below: 


T o = T D * (w d /w) 
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where = 6. 25 seconds is the value of r Q at the design flow rate. An equally impor- 
tant result was found that the lag function constant L was given by the relation 

L = l d * ( w d /w > 


where 


L D = 50. 0 sec 

These findings make the results more general; they can be extended to specific flow 
rates not treated herein. Very good agreement is shown in figure 2 between the Bode 
plots of the simplified transfer function and the computer program for inlet temperature 
disturbance over the frequency range studied. 

The dynamic response of the SNAP- 8 radiator exit temperature to a change in inlet 
flow rate is indicated by the Bode plots shown in figure 3. The assorted symbols show 
the computer program generated points. These points display a small resonance effect 
in the midfrequency range; however, to account for this effect would overly complicate 
the transfer functions. Asymptotes fitted to the data points indicated that the simplified 
transfer function should be a first-order lag. The phase angle approaching 90° at higher 
frequencies also points to a first-order lag form of transfer function. Therefore, the 
simplified transfer function for outlet fluid temperature to fluid flow rate was chosen to 
be 


(ifl°w) ln (! + T ^S)1 

Flow rates of 10, 20, 50, and 100 percent of design point flow rate are shown in fig- 
ure 3(a) to (d), respectively. Again, the important relation for the time constant shown 
below is valid 


t 4 = t D* (w d /w > 


where 


Tp = 45. 0 sec 

is the value of at the design flow rate. The solid lines in figure 3 show the locus of 
points given by the simplified transfer function model. 
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The value of time constant shown was chosen to give a good average fit over the fre- 
quency range studied. The agreement between the simplified transfer function and com- 
puter generated results can be improved for lower frequencies by varying the time con- 
stants for specific cases where this may be required. 

Since the simplified transfer functions are given in normalized form, they must be 
combined with a steady- state gain factor in their application. The steady- state ampli- 
tude ratios of fluid outlet temperature to inlet fluid temperature and to inlet flow rate 
are plotted against flow rate in figure 4. 

SUMMARY OF RESULTS 

From the analysis of the dynamic characteristics of a typical waste-heat radiator 
for SNAP- 8, the following results were obtained: 

1. For the assumed radiator, the transfer function relating changes in outlet tem- 
perature to changes in inlet temperature is close to a third-order lag multiplied by a 
dead-time element. The time constant of the third-order lag varies inversely with flow 
rate. The design flow rate time constant is 6. 25 seconds. The dead time also varies 
with flow rate, being equal to the design point value (50 sec) multiplied by the ratio of 
the design point flow rate to the given flow rate. This type of transfer function agreed 
very well with the frequency response generated by the computer simulation of the radi- 
ator for flow rates from 10 to 100 percent of design flow. 

2. The simplified transfer function for the radiator relating changes in outlet tem- 
perature to changes in flow rate is a first-order lag. The time constant of the first- 
order lag varies with flow rate, being equal to the design point time constant (45. 0 sec) 
multiplied by the ratio of design point flow rate to the given flow rate. A resonance ef- 
fect in the frequency response generated by the computer simulation precluded the sim- 
plified transfer function from matching this frequency response closely in the midfre- 
quency ranges. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 13, 1968, 

701-04-00-02-22. 
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APPENDIX - COMPUTER PROGRAM LISTING 


$ I B FTC MAIN P LIST 

THIS IS A FORTRAN IV SOURCE DECK PROGRAM FOR 
COMPUTING THE FREQUENCY RESPONSE OF THE SNAP-8 
RADIATOR FOR PERTURBATIONS IN INLET TEMPERATURE 
OR INLET FLOW RATE. THE ASSUMED RADIATOR DESIGN 
IS A CYLINDER 9.260 METERS IN LENGTH AND 3.72 
METERS IN DIAMETER. THE PROGRAM IS GENERAL 
EXCEPT FOR THE INITIALIZATION OF CONSTANTS. 

METRIC UNITS ARE USED. 

COMMON TIME,DTIME,C7,FL0WD,TDRA(41 ) ,TRA(41 ) , F I NEF , DF I N1 A , THE AT 

700 FORMAT ( 1H 12F10.5) 

701 FORMAT (94H1 TIME TDRA ( 41 ) TRA(41) FINEF DF 

1IN14 TRAD TIMIL THEAT ) 

702 FORMAT (1H ,8F12.5) 

T I ME=0 .0 
CT I ME=0.0 
DT I ME=0 • 1 
C 7 = 6 1 9 . 

FLOWD=. 454/153.0 
WRITE (6.701 ) 

10 CALL RADTOR 

CTIME=CTIME+DTIME 
I F ( C T I ME-4 . 9 5 ) 500,30,30 
30 CTIME=CTIME-5.0 
T=0 . 0 

DO 35 1=2,40 
35 T=T+TRA ( I ) 

TRAD=( (TRA(l)+TRA(4l) ) /2 . 0+T ) /40 . 0 

WRITE (6,702 ) T I ME , TDRA ( 41 ) , TRA ( 41 ) , F I NEF , DF I N 14, TRAD, T I M I L ,THE AT 
500 TIME=TIME+DTIME 

IF(TIME— 1200.0) 10,10,3 

3 TMAX=TDRA ( 41 ) 

TM I N=TMAX 
TIME=. 00001 
DT I ME=0 • 1 1 

READ (5,700) AMPL , OMEGA, STEP 
WRITE (6,700) AMPL , OMEGA, STEP 
12 I F ( T I ME • LE • 50 • ) AMPLU = AMPL*T I ME/50 . 

T DR A ( 1 ) = 619. +AMPLU*SIN(OMEGA*TIME)+STEP (Eq. (13)) 

C OR FL0WD=(4.54+AMPLU*SIN(0MEGA*TIME)+STEP)/153.0 (Eq. (14)) 

CALL RADTOR 
CTIME=CTIME+DTIME 
IF(CTIME — .95) 504,34,34 

34 CTIME=CTIME— 1.0 

WRITE (6,700) TIME, THEAT, FLOWD, (TDRA( I ), 1=1, 41, 5) 

504 T I ME=T I ME+DTIME 

T MAX = AMAX 1 ( TMAX , TDR A ( 41 ) ) 

TM I N=AM I N1 ( TM IN, TDRA (41 ) ) 

IF(TIME— 2500.0) 12,12,5 

5 AMPLR=(TMAX-TMIN)/(2.0*AMPL) 

WRITE (6,700) AMPLR 

STOP 

END 
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$ I BFTC RADI AT LIST 

SUBROUTINE RADTOR 

COMMON TIME,DTIME,C7,FL0WD,TDRA(41 ) ,TRA(41 ) , F I NEF , DF I N14 , THE AT 
DIMENSION HEAT (41 ) ,FINF(22 ) ,C0N(2 ) 

DATA ( F I NF ( I ) » I =1 , 2 2 ) / 1 . 0, . 854 , . 78 3 , . 72 7 , . 68 1 , . 642 , . 608 , . 578 , . 554 , 

1.532, .512, .495, .480, .466, .453, .441 , .430, .420, .410, .401, .393, .385/ 

DATA A159,RETOT , CON ( 1 ) ,CON( 2 ) , AARM , AF I N , C 1 , F I NW I , F I NTH , 

1 HCAPTA ,HCAPF/ 1 70 . , . 02 1 84 , 2 33 .3 , . 0521 2 , . 00 1 9 1 8 , .0 1 5 , 5 . 103E-8 , 

2 .03414, .000762,35.26, 1.977/ 

THEAT=0 .0 

IF(TIME)502,502,503 

502 FAC TOR = ( 1 »-EXP( —CON ( 2 ) / FLOWD ) ) / ( CON ( 2 ) / FLOWD ) 

INITIAL CONDITIONS 
DO 510 1=1,41 
TRA ( I ) =C 7 
510 T DRA ( I ) =C 7 
RETURN 

MAIN BODY OF SUBROUTINE 

503 DO 520 J = 1 ,41 

SI=C1*TRA( J )**3*FINWI**2/ ( A159*FINTH) 

I=SI*10. 

IF(I.LT.O) 1=0 
S I=S I-FLOAT ( I > / 1 0 . 

B=( F INF U + l ) -F INF < 1+2) 1*10.0 
F I NEF=F I NF ( 1+1 >-B*SI 
FIN14=S0RT ( SORT ( FINEF ) ) 

DFIN14=-B*3.0*S I / ( F I N14**3*4 . 0*TR A ( J ) ) 

INTEGRATION 

HEAT (J )=C1*(AARM+AFIN*FINEF)*TRA(J)**4 

TRA ( J ) =TRA ( J ) + ( ( TDRA ( J ) -TRA ( J ) ) *F ACTOR/RE TOT-HEAT ( J ) ) / ( HC APTA + HC AP (Eq. (12)' 
1F*( FIN14+TRA ( J )*DFIN14) ) *DT IME 
520 THEAT=THEAT+HEAT (J) 

DXR = DXR + CON ( 1 ) *FLOWD*DT I ME 
I F ( DXR- 1 . >505,550,550 
550 DXR=DXR-1. 

FACT1=1.— EXP(— C0N(2 ) /FLOWD ) 

FACT0R=FACT1/ (CON (2 ) /FLOWD) 

DO 560 1=1,40 
K=42- I 

560 TDRA(K)=TDRA(K-1 )+FACTl*( TRA(K— 1 )-TDRA(K-l ) ) (Eq. (5)) 

IF ( DXR-1 . ) 505,550,550 
505 RETURN 
END 
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Figure 1. - Assumed design of SNAP-8 radiator. 
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Figure 2. - Frequency response of fluid outlet 
temperature to fluid inlet temperature for 
various percentages of design flow rate. 
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D Computer program 



(a) Flow rate of 10 percent of design value. 




(D) Flow rate of 20 percent of design value. 
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(c) Flow rate of 50 percent of design value. (d) Flow rate of 100 percent of design value. 

Figure 3. - Frequency response of fluid outlet temperature to inlet flow rate for various percentages of 
design flow rate. 
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(a) For perturbations in inlet temperature. 
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(b) For perturbations in inlet flow rate. 

Figure 4. - Steady-state amplitude ratios against flow rate 
for perturbations in inlet temperature and inlet flow 
rate. Design flow rate, 4.54 kilograms per second; 
inlet fluid temperature, 619.0 k. 
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